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Abstract 

We analyze the thermal conductivity of frustrated spin-1/2 chains within linear-response theory focusing on its 
frequency dependence at finite temperatures. Using exact diagonalization, the low-frequency limit of the thermal 
conductivity is discussed in the high-temperature regime. 
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The recent interest in thermal transport properties 
of one-dimensional spin-1/2 systems is motivated by 
the experimental observation of significant contribu- 
tions to the thermal conductivity of quasi-one dimen- 
sional materials originating from magnetic excitations 
(see, e.g., Refs. [1,2] and further examples mentioned 
in Ref. [3]). A number of recent theoretical papers have 
either focused on the integrable XXZ model [3,4,5,6,7] 
or on generic spin models like frustrated chains or spin 
ladders [3,6,8,9,10,11,12]. When analyzing transport 
properties, the real part of the conductivity is usually 
decomposed into a singular part at zero frequency uj = 
with weight Dth and a regular part 



ReK{uj) = Dth{T)3{uj) + Krceioj) . 



(1) 



Here, k denotes the thermal conductivity, Dth is the 
thermal Drude weight, T is the temperature, and S{u)) 
is the (5-function. 

The XXZ model shows ballistic thermal transport 
properties due to the exact conservation of the energy- 
current operator [4] . A complete understanding of the 
temperature and magnetic field dependence of the rel- 
evant quantity, the thermal Drude weight Dth, has re- 
cently been established [3,5,6,7]. 
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Regarding spin ladders and frustrated chains, both 
numerical [3,6,12] and field theoretical studies [3,11] 
indicate normal transport behavior, characterized by a 
vanishing thermal Drude weight in the thermodynamic 
limit. Thus, the information relevant for the interpre- 
tation of experiments is encoded in the low-frequency 
behavior of the thermal conductivity. 

This paper focuses on the example of the frustrated 
chain in the massive regime (see Ref. [12] for a similar 
study of spin ladders) . The Hamiltonian reads 



ii' = ^ ft; = J^[S, • 5^+1 -f a5, • Si 



l+2\ 



(2) 



hi is the local energy density; Si is a spin-1/2 operator 
acting on site I. N denotes the number of sites and 
periodic boundary conditions are imposed. aJ is the 
next-nearest neighbor coupling; for a = 0, we obtain 
the spin-1/2 Heisenberg chain. In the following, we will 
consider the case of a = 1 . The energy-current operator 
corresponding to this Hamiltonian is chosen such that 
the equation of continuity —dthi = divjth i is fulfilled 
[6] 

N N 1 

jth = J^ jth,! = «X1 XI [hl-r-l,hl + a] . (3) 



i=l 



! = 1 r,s=0 



Following linear-response theory, the thermal conduc- 
tivity can be written as [13] 
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Fig. 1. Thermal conductivity of frustrated chains with a — 1 
[see Eq. (1)]. Panel (a): regular part of Re k{uj) vs uj for 
T/J ^ 1, 2 for TV ^ 18 and N ^ 20 sites. Panel (b): integrated 
weight I{uj) [see Eq. (5)] vs u for various temperatures. Inset 
of (b): T ■ Jo vs temperature T. 
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J2 Pn\{m\j,i,\n)f5{ij-AE) (4) 



with p„ = exp(— _E„/T)/Z being the Boltzmann 
weight, Z the partition function, AE = Em — E^, and 



ai^) = [1 



,-^/T- 



/uj. En are eigen-energies and 



eigen-states of the Hamihonian Eq. (2). Another useful 
quantity is the integral of Re k{ijj) over frequency u) 

liuj) = Dti,{T) + 2 [ duj' K,eg(c^'); ^0 = hm I{uj) .(5) 



0+ 

Applying standard numerical techniques allows us to 
study systems with A*' < 20 sites. The dimension of the 
largest subspace is 9250 after exploiting translational 
invariance and conservation of total S^ . 
Let us now turn to the discussion of the results pre- 
sented in Fig. 1. Panel (a) contains the regular part 
Krcg(tj) for T/J = 1, 2 and A = 18 and A = 20. In this 
plot, Krog(<^) has been binned with AE/J ^ 10~^. 

First, note that n,-cg{u;) is a structureless function of 
frequency u;. Depending on temperature, the data are 
well converged down to frequencies to/ J < 0.25 while 
deviations between A = 18 and A = 20 sites become 
visible for lower frequencies. Furthermore, we observe a 
characteristic down-turn of Krcg(<^) at low frequencies 
which is also present for spin ladders [12]. This feature 
is less pronounced on larger systems because spectral 
weight is transferred from the Drude weight to finite 
frequencies as the system size increases. Moreover, we 
have to keep in mind that the Drude weight is excluded 
in the curves displayed in Fig. 1. Note that the Drude 
weight amounts to less than 3% of the total weight Iq 



from Eq. (5) for A = 20 sites and T/J = 1. Thus, the 
scenario of ballistic thermal transport is unlikely for 
frustrated chains. We refer the reader to Refs. [3,6] for a 
detailed scaling analysis of the Drude weight. While 
the down-turn at low frequencies renders it difficult to 
draw precise conclusions about the functional form of 
the frequency dependence, we can, however, fairly well 
estimate the absolute value of k = lim„-^o Krcg(<^). For 
instance, the data shown in Fig. 1 (a) indicates k/ J^ ~ 
(0.27 ± 0.3) for T/J = 2. 

Next we comment on the integrated spectral weight 
I{uj), Eq. (5), which is shown as I{u))/I(, in Fig. 1 (b) 
for A = 20 sites and various temperatures T/J > 0.3. 
Here, we see that first, I(lu)/I(, approaches 1 the faster 
the larger T is. Second, the curve for T/J = 0.3 is less 
smooth than the others indicating the limitations of 
an exact diagonalization study with respect to the ac- 
cessible temperature range. A smooth curve ReK(aj) is 
obtained for T/J > 0.5 for A = 20. Finally, we men- 
tion that 7o ~ T~^ for high temperatures, as expected 
for thermal transport [see inset of Fig. 1 (b)]. 
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